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The nuclear wobbling motion in the Lu region is studied by the microscopic cranked 
mean-field plus RPA method. The Woods-Saxon potential is used as a mean-field with 
a new parameterization which gives reliable description of rapidly rotating nuclei. The 
prescription of symmetry-preserving residual interaction makes the calculation of the 
RPA step parameter-free, and we find the wobbling-like RPA solution if the triaxial 
deformation of the mean-field is suitably chosen. It is shown that the calculated out- 
of-band B(E2) of the wobbling-like solution depends on the triaxial deformation in the 
same way as in the macroscopic rotor model, and can be used to probe the triaxiality 
of the nuclear mean-field. 



1. Introduction 

In this talk we would like to discuss a fundamental question: How nucleus rotates if 
it is triaxially deformed? Namely, the question is whether nucleus exhibits the so- 
called wobbling rotational motion. We have been conducting research on this subject 
for many years, and we would like to report, especially, on the recent development 
of our microscopic study of the wobbling rotation. 

The wobbling motion is an analogy of the quantized motion of the macroscopic 
("rigid-body") triaxial rotor. The spectra of the model hamiltonian of triaxial rotor 
can be solved approximately (using the H — 1 unit) like this; 1 

I 2 I 2 I 2 1(1 + 11 

H mt = ^ + ^ + ^, => B (j,„)«£i+J + £( ^, b („ + i/2), (1) 

Ux Uy U z £\J x 

where the wobbling energy w wo b is given by 



W wo b — 




Here it is assumed that the main rotation axis is the a;-axis {J x > J V tJz). This 
type of rotational motion is only possible when the system is triaxially deformed, 
and appears as a multiple band structure shown in Fig. 1. It is composed of a 
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group of bands lying on top of each other, whose horizontal sequences are usual 
rotational bands connected by the strong stretched (\AI\ = 2) E2 transitions, while 
the vertical sequences are connected by slightly weaker \AI\ — 1 E2 transitions. In 
this way, the vertical excitation, which corresponds to tilting the angular momentum 
vector (in the body-fixed frame), can be regarded as the phonon-like excitation, 
n = 1, 2, ... in Eq. (1), whose phonon energy is given by the famous formula (2). 



R 



Wobbling 




3-phonon 
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Fig. 1. A schematic illustration of the band structure of the wobbling motion. 



Until quite recently, there is no definite evidence of the multiple band structure 
that exhibits this characteristic feature. It has, however, been first discovered in 
163 Lu: 2 The one phonon wobbling band has been measured nowadays in some Lu 
isotopes, 3 ' 4,5 and moreover the two phonon wobbling band has been identified in 
163 Lu. 6 These bands are examples of the so-called triaxial superdeformed (TSD) 
bands in the Lu and Hf region. The measured in-band (AI = —2) B(E2) values 
between the horizontal sequences in these TSD bands are typically about 500-700 
Weisskopf units, and the out-of-band (AI = —1) B(E2) values from the one-phonon 
wobbling band to the yrast TSD band are about 100 units or more. 7 These values 
are very large and consistent to the wobbling picture predicted by the macroscopic 
particle- rotor model, 8,9 where an odd 213/2 quasiproton is coupled to the rotor. 

The triaxial deformation is crucial to realize the wobbling motion. Then, it is 
important to know how the effects of triaxiality appear in observables. One is the 
B(E2), which are related to the two intrinsic quadrupole moments, Q20 and Q22, 



Q 



20 




]T(2z 2 - x 2 - y 2 ) a , Q 22 



15 
327 



and the triaxiality parameter 7 is defined in terms of them as usual: 

\/2(Q 22 ) 



tan 7 



(Q 



20/ 



(Lund convention of sign). 



(3) 



(4) 



The other is the energy spectra, which reflects the properties of the three different 
moments of inertia about intrinsic axes. Here, it should be pointed out that the 
characteristic feature of the out-of-band E2 transitions observed in Lu nuclei sug- 
gests the "positive 7" shape, i.e. the main rotation axis (the x-axis in our notation) 
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is the shortest axis in that shape. This shape completely contradicts the well-known 
irrotational moments of inertia, in which J y is the largest, and then the wobbling 
frequency (2) becomes imaginary. Therefore, we don't know what kind of inertia 
should be used in macroscopic models, and a microscopic approach is necessary. 

2. Microscopic approach 

In order to understand the wobbling motion in Lu nuclei, we use a standard mi- 
croscopic approach; the cranked mean-field and the random phase approximation 
(RPA). The basic idea is to describe the horizontal rotational bands in Fig. 1 by the 
semiclassical cranking prescription, and at the same time, the vertical phonon-like 
excitation in terms of the RPA, which is known to be successful for describing the 
collective vibrational modes. Namely, 

ti = h dcl - uJrotJx, (cranking) (5) 
[H',Xl]=u n Xl (RPA) withtf' = /i' + Mnt, (6) 

where hdef is the triaxially deformed mean-field hamiltonian with the cranking 
frequency w ro t about the x-axis, X^ is the creation operator of the n-th RPA eigen 
mode with the eigen energy u) n , and Vi n t is the residual interaction. 

It should be emphasized that the residual interaction is constructed so as to 
restore the rotational symmetry broken by a general mean-field, i.e. 

[hdef, Jk}^Q => [hdef + V int} Jk] = 0, (k = x, y, z), (7) 

with 

VSnt = -^ ^2 KkF k> F k = [h dcf ,iJ k ], K k = ([[hdef,Jk],Jk\), (8) 

k— x,y,z 

where the expectation value is taken with respect to the cranked mean-field yrast 
state. Thus the RPA solutions are determined without any ambiguities by a given 
mean-field, and there are no adjustable parameters in the RPA step. 

However, it is not clear how the result of such a microscopic approach is re- 
lated to the wobbling picture of the macroscopic rotor model. It was shown by 
Marshalck 10 many years ago that, if an appropriate RPA mode exists, the rotor 
model wobbling picture naturally appears by going from the uniformly rotating 
(UR) frame, where the cranking prescription is applied, over to the principal axis 
(PA) frame. Taking into account the fact that the RPA can be regarded as a small 
amplitude limit of the time-dependent mean-field theory, the time-dependent mean- 
field describing the wobbling motion is given in the UR frame as 

hvR(t) = hdef - u rot J x - K, y Ty{t) Fy - k 7 T ' z (t) F z , (9) 

where Tk{t) = (-Ffc)uR.(i) (k — y 7 z) are the UR frame expectation values; only 
the y and z components are relevant to the wobbling phonon excitation, which 
transfers spin by AI = ±1. In the case of the harmonic oscillator potential, e.g. 
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the Nilsson potential, the operators F y and F z are proportional to Q y and Q z , the 
non-diagonal quadrupole tensor operators. Even in the case of general potentials, 
these non-diagonal parts Qk {k = x, y, z) are used to define the PA frame; 

A15 A 

(Qfe)pA = 0, (PA condition), Q k = J — y ^{ x i x j)a, (k,i, j)-cyclic. (10) 

a=l 

The shape fluctuation vanishes by this PA frame condition, but the fluctuation 
of the angular momentum vector appears in place of it; thus the time-dependent 
mean-field in the PA frame is transformed as 

hpA(t) = h dci - LU x (t)J x - LUy(t)J y - U) z (t)J z , (11) 

where uj x (t) ss aj rot in the small amplitude limit, and the presence of u> y (t) and 
u> z (t) reflects that the angular frequency vector also fluctuates about the main 
rotation axis (the z-axis) . The relation between the UR and PA frames is depicted 
schematically in Fig. 2. The PA frame corresponds to the body- fixed frame of the 
rotor, and the moments of inertia are naturally introduced by 

Jx = (Jx)/w rot , J y {n) = Jy{n)/uj y {n), J z (n) = J z (n)/u> z (n), (12) 

where Jfe(n) and LOk(n) (k = y, z) are Fourier components of the PA frame expecta- 
tion values with respect to the n-th RPA cigen mode. Using these microscopically 
calculated moments of inertia, the RPA excitation energy can be written in exactly 
the same form (2) as in the macroscopic rotor model. It can be also shown 11 that 
the out-of-band B{E2) calculated by the RPA approach can be expressed in the 
same way as in the rotor model. 




y y 



Fig. 2. A schematic illustration depicting the relation between two dynamical pictures in the 
uniformly rotating (UR) frame and the principal axis (PA) frame. 



3. Results of microscopic calculations 

Realistic calculations based on the Marshalek's theory have been carried out in the 
previous works. 12 ' 11 After the discovery of the wobbling phonon excitation in the 
TSD bands in the Lu and Hf region, new calculations have been performed, and 
it is confirmed that the wobbling-like RPA solutions do exist. 13 ' 14 However, the 



9, 2008 0:42 WSPC/INSTRUCTION FILE shanghai06 



Instructions for Typing Manuscripts (Paper's Title) 5 

Nilsson potential has been used in these calculations, and they are suffered from 
the problem of large moments of inertia because of the spurious velocity dependence 
coming from the Z 2 -term. Therefore, we have conducted new research with using the 
Woods-Saxon potential as a nuclear mean-field. Especially, a new parameterization 
of the Woods-Saxon potential has been provided quite recently by Ramon Wyss, 
with which it is possible to nicely reproduce the basic properties like the neutron 
and proton density distributions. This is very important to obtain reliable results 
for the moments of inertia and the quadrupole moments. 

3.1. Mean-field parameters 

As it is discussed in the previous section, there is no adjustable parameter in the 
RPA step once the mean-field is specified. Therefore, let us briefly explain the mean- 
field parameters used in our calculation. In principle, they should be determined 
sclfconsistently by minimizing the energy. Actually we are developing the Woods- 
Saxon Strutinsky calculations, but it is not available yet. Therefore we have taken 
the deformation parameters corresponding to the minimum of the Nilsson Strutin- 
sky calculations, an example of which is shown in Fig. 3. The TSD minimum moves 
slightly as a function of the rotational frequency but the amount of change is small, 
so that we have used average values, /?2 = 0.42, 7 = 12°, /3 4 = 0.034 in all the 
calculations, we have checked that the qualitative feature of the calculated results 
does not change if these values are modified in a reasonable range. 



10 10 
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Fig. 3. An example of the potential energy surface obtained by the cranked Nilsson Strutinsky 
calculation for the 7ri 13 /2 (k,ci) = (+,+1/2) configuration at / = 53/2+ in 163 Lu. The energy 
between contours is 250 keV. In the left figure is adopted the "usual" 7, which is used to parame- 
terize the Nilsson potential, while in the right one the calculations are same but the 7 defined by 
a similar Eq. to (4) (see the footnote in the text) is adopted. 

Here we would like to point out that there are various different definitions for 
the triaxiality parameter 7. In this presentation, we use the one defined by Eq. (4) 
throughout, which is directly related to the two intrinsic quadrupole moments, and 
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so the B(E2). In Fig. 3, the same calculations are employed but two representations 
with different definitions of 7 are shown. The left part is an usual representation, 
where the 7 is the one used in the Nilsson potential, and the TSD minimum has 
7 w 20° and 62 w 0.4. In the right part the definition similar to Eq. (4) is used a keep- 
ing the £2 as usual, and then the TSD minimum has 7 w 12°. Thus, the difference 
between the two definitions for the same shape is rather large for large €2 deforma- 
tion like in the case of the TSD bands. More details about the various definitions 
of 7 parameters and relations between them will be discussed elsewhere. 15 

The neutron and proton pairing gaps are also important mean-field parameters. 
We have used the following parameterization, which is convenient to avoid the 
abrupt pairing collapse; 



Z\(w rot ) = A) x 



2\uj c 

2 V^rot 



^rot < ^cj 



^rot > L^c, 



(13) 



where Ao is given by the even-odd mass differences of neighboring even-even nuclei, 
and uo c is determined by the selfconsistent monopole pairing calculation. The resul- 
tant pairing gap parameters are shown in Fig. 4. Now all the mean-field parameters 
are fixed, and there is no adjustable parameter for the following RPA calculations. 




3.2. RPA calculations 

Now let us present the results of RPA calculations for the TSD wobbling band in 
163 Lu, for which most extensive data are available. 7 ' 16 ' 17 First, Fig. 5 shows the 
cranking moment of inertia, J x in Eq. (12), as a function of the rotation frequency. 

a More precisely, the expectation values in Eq. (4) is replaced by those with respect to the liquid- 
drop like sharp cut-off density distribution. Then the 7 is determined purely by the geometry of 
the potential without recourse to actual wave functions. 
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This inertia is related to the slope of the horizontal rotational band in Fig. 1. Our 
calculation agrees the trend of experimental data rather well. As for a reference, 
here and in the following, we also include examples of the particle-rotor model 
calculations by Hamamoto-Hagemann; 9 the relevant parameters of the model are 
7 = 20°, and J~j R) , j£ R) , J~j R) = 48,45,17 T^MeV" 1 , respectively. The result of 
the particle-rotor model is monotonically decreasing, which can be understood by 
the following approximation valid for a highly aligned quasiparticle band; J x 
Jx + j'/^rot, where j ~ 13/2 is the maximal alignment of the 7rz 13 / 2 orbit. 



100 




40 1 ■ 1 ' 1 ■ 1 ■ 1 ' 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 
hco rot [MeV] 

Fig. 5. The moment of inertia, J x in Eq. (12) as a function of the rotational frequency in 163 Lu. 
The result of the particle-rotor model is also included as the dotted line with the legend "p-rotor" . 
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Fig. 6. The excitation energy of the one-phonon wobbling band as a function of the rotational 
frequency in 163 Lu. 



Next we show the results of the RPA. We have found that the wobbling-like RPA 
solution does exist also for the Woods-Saxon mean-field potential; the calculated 
excitation energy is compared with experimental data in Fig. 6. The solution is 
stable against the change of mean-field parameters in a reasonable range, so we be- 
lieve that the existence of the wobbling-like solution has been confirmed. Although 
the result qualitatively agrees, the calculated excitation energy is smaller than the 
measured one, and it vanishes at about w rot = 0.52 MeV. For the excitation energy, 
the particle-rotor model gives increasing energy as a function uj Iot in contrast to the 
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data. It is, however, difficult to change this general trend as long as the constant 
moments of inertia of the rotor are used, see Eq. (2). 
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Fig. 7. There moments of inertia calculated by the Marshalck's theory, Eq. (12), applied to the 
wobbling-like RPA solution in 163 Lu. 



In Fig. 7, we show three moments of inertia calculated according to the Mar- 
shalck's theory by Eq. (12). Note that they are not the inertia of the rotor, but 
the total inertia of the system; it is not possible to divide the contribution into 
the particle and the rotor parts a priori in the microscopic calculation. They are 
not constant but gradually change, which is necessary to understand the decreasing 
trend of the observed excitation energy 





Fig. 8. The out-of-band over in-band B(E2) ratio (left) and the in-band B(E2) (right) as a function 
of the rotational frequency for the one-phonon wobbling band in 163 Lu. The lines connected with 
star symbols are results of the specific RPA calculation, where the mean-field parameter 7 is 
changed from 7 = 12° at ui Iot = 0.2 MeV to 7 = 22° at w rot = 0.4 MeV. 



Now we come to one of the most important observables, the B(E2), which are 
depicted in Fig. 8. The measured out-of-band over in-band B(E2) ratio is almost 
constant (the left part of Fig. 8), but a gradual decrease of the B(E2) ratio as 
a function of u> rot is expected in the rotor model if the triaxiality parameter 7 is 
kept constant. Both ratios calculated by the RPA approach and by the particle- rotor 
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model follow this tendency, but our microscopic result is about 60% of the measured 
value at iv Iot w 0.3 MeV, where the particle-rotor model gives correct magnitude. 
However, the calculated out-of-band B(E2) values are more than 50 Weisskopf 
units, which is a huge value for a normal RPA solution, and clearly indicates that 
the obtained solution is not of vibrational but of rotational character (B(E2) values 
of typical collective vibrations in deformed nuclei are about 10 Weisskopf units at 
most). The average values of the B(E2) ratio in the particle- rotor model agrees 
rather well; in fact the parameter 7 = 20° is chosen for this reason. 9 Note that our 
mean-field value is 7 — 12°, and is much smaller; we have checked that almost the 
same larger value as in the particle-rotor model can be obtained if we use a similar 7 
value. As for the in-band B(E2), it is just related to the quadrupole moment of the 
mean- field about the rotating axis (rc-axis). The calculated result (the right part of 
Fig. 8) is almost constant because we have used constant deformation parameters, 
while the data show a clear tendency to decrease. 

As is discussed already, the measured B{E2) strongly suggest that the deforma- 
tion of the mean-field, especially the triaxiality 7, is changing as a function of the 
rotational frequency or spin. In the simple rotor model, both the in-band B(E2) 
and the B(E2) ratio can be estimated approximately as 

B(E2) in » ^|e 2 £> 2 - z 2 )) 2 = J-e 2 Q 2 cos 2 ( 7 + 30°), (14) 



B(E2) out ^ 2 
B(E2) in ~ / 

with 



(w z /w y y/ 4 sin(7 + 60°) + (w y /w z ) 1 / 4 sin 7 
cos(7 + 30°) 



(15) 



Wy = 1/J Z ~ 1/J X , W Z = l/Jy - 1/J X . (16) 

The in-band B(E2) depends on both the magnitude of deformation, e.g. e 2 , and 
the triaxiality 7, while the B(E2) ratio depends only on 7 in this approximation. 
Moreover, these expressions clearly shows that the B(E2) are quite sensitive to the 
triaxiality parameter 7. The 1/1 factor in Eq. (15) gives the decreasing trend of 
the B(E2) ratio if 7 and J x ,J y ,J z are kept constant. Now, can we understand 
the measured spin dependence of B(E2) by changing the 7 parameter? We have 
tried to play a game by increasing 7 linearly from 7 = 12° at co TOt = 0.2 MeV to 
7 = 22° at oj Iot = 0.4 MeV, with keeping other parameters unchanged. The results 
are shown by the lines with star symbols in Fig. 8. Both the B(E2) ratio and the 
in-band B(E2) nicely agree with data. 



4. Summary 

We have performed microscopic RPA calculations using the Woods-Saxon potential 
as a mean-field. We have found that the wobbling-like eigen mode does exist in 
this calculation. Considering that the similar solutions have been obtained in the 
previous investigation using the Nilsson potential, 13 ' 14 we believe that this confirms 
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the existence of the collective mode, the nuclear wobbling motion, suggested by the 
macroscopic rotor model. 

The calculated excitation energy is smaller than the experimental one by about 
100-150 keV. This may suggest that some part of coupling effects to the odd «i3/2 
quasiproton is missing in our calculation. In this respect, Hamamoto-Hagemann 
has discussed, 9 by using a simple approximation in the particle-rotor model, that 
the effect of coupling to the odd particle appears as a constant energy shift to 
the wobbling energy (2), Aj wo b = j/Jx, where j is the quasiparticle alignment. 
Applying this estimate to the case of 163 Lu, j = 13/2 and J x w 70 fi 2 MeV _1 , the 
shift amounts to Au wo b ~ 93 keV. This fills most part of the energy underestimated 
in our calculation, although we are not so sure whether such an energy shift can 
be justified in our microscopic framework: An explicit particle-vibration coupling 
treatment based on our framework is necessary in order to elucidate this point. 

The experimental out-of-band over in-band B(E2) ratio on average suggests 
larger triaxiality 7 s=s 20° than the value obtained for TSD minima in the Nilsson 
Strutinsky calculation, 7 w 12°; there is an apparent discrepancy. Note that 7 
here is defined by Eq. (4), and not the "usual" 7 in the Nilsson potential, whose 
translated value is incidentally about 20° in the case of the TSD bands. 

Both the in-band B(E2) and the B(E2) ratio indicate an increase of the triaxial 
deformation as a function of the rotational frequency or spin. We can obtain nice 
agreements between the microscopic RPA calculations and the experimental data, 
if the triaxiality parameter 7 is changed linearly, e.g. from 7 = 12° to 22°. However, 
the Nilsson-Strutinsky calculation suggests almost constant 7 values for TSD min- 
ima, and it may be difficult to expect such a large change of 7 values. We definitely 
need more study for quantitative understanding of the observed properties of the 
nuclear wobbling motion in the Lu region. 
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